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INTRODUCTION 


The  circuit  and  systems  analysis  program  SUPER*SCEPTRE  is  a  powerful  and 
highly  Interdisciplinary  simulation  language.  Among  its  various  attributes  is 
Its  ability  to  utilize  user  defined  models  as  circuit  or  system  elements.  Unis 
report  describes  a  procedure  for  automating  the  creation  of  diode  and  bipolar 
transistor  (NPN)  models  used  by  SUPER*SCEPTRE.  While  the  final  models  generated 
appear  in  SUPER*SCEPTRE  format,  they  can  he  easily  adapted  to  the  format  of  any 
nonlinear  circuit  simulation  program.  The  modeling  procedure  mentioned  begins 
with  the  acquisition  of  semiconductor  data  by  means  of  an  automated  test  setup. 
The  setup  Is  controlled  by  a  Hewlett  Packard  9816S  microcomputer  in  concert  with 
an  H.P.  4145B  semiconductor  parameter  analyzer.  The  data  is  then  transformed  by 
the  9816S  Into  a  SUPER*SCEPTRE  model.  A  discussion  of  the  parameters  used  to 
describe  each  model  will  first  be  presented  followed  by  a  description  of  the  test 
configurations  used  to  collect  parameter  data. 

DIODE  MODEL 

The  diode  model  for  which  the  following  parameters  will  hi?  developed  is 
taken  from  reference  l  and  appears  in  figure  l.  Diode  model  parameter  defini¬ 
tions  are: 

CD  ■  The  sian  of  the  diode  transition  and  diffusion  capacitances  whore: 


and 


Transition  capacitance 


<phi-vcd) 


Co 

—  N  “  01 


Diffusion  capacitance  «  KD  (.JD  +  !<.) 


JD  -  Current  generator  representing  the  diode  junction  current, 
generator  is  a  function  of  the  voltage  VC1). 

JD  -  Is  [EXP( Theta  VCD)  -  lj 

I,,  *  Diode  reverse  saturation  current 

Theta  -  Constant  of  the  diode  equation 

RB  -  Diode  bulk  resistance 

Co  »  Constant  of  the  transition  capacitance  equation 
-  Junction  contact  potential 
■  Junction  grading  constant 


This 


PHI 

N 

KD 
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_ , , ,  .  .  _  _  Theta 

■  Diffusion  capacitance  constant  *  r- — — —  — r- 

(  6  •  2  o  F) 

-  Frequency  parameter 


Vqj  ■  Voltage  across  capacitor  CD,  which  Is  equal  to  the  diode  junction 
voltage 

Diode  Equation  Constant  and  Reverse  Saturation  Current 

The  diode  equation  constant  Theta  is  defined  as 

Th*t*  ■TkV  (1) 

where  q  is  the  charge  of  an  electron,  k  is  Botzman's  constant,  and  T  Is  the 
diode's  junction  temperature  in  Kelvins.  Theta  is  obtained  from  the  diode  junc¬ 
tion's  measured  forward  biased  characteristic  and  the  assumption  that  the  char¬ 
acteristic  can  be  closely  approximated  by 

JD  *  Is  EXP  (Theta  VCD)  (2) 

where  I<,  represents  the  diode's  saturation  current  and  VCD  1  ts  forward  junction 

voltage.  Diode  D' s  forward  characteristic  Is  measured  by  using  figure  2.  The 
forward  Junction  current  and  the  foward  junction  voltage  are  recorded  by  means  of 
source/monitor  units  SMU3  and  SMJ1,  respectively,  for  automatically  varied  cur¬ 
rent  values  of  SMU1.  Source/monitor  units  are  an  integral  part  of  the  HP4145B 
and  act  as  either  voltage  sources  and  current  monitors  or  current  sources  and 
voltage  monitors.  In  this  test,  SMUl  consists  of  a  current  source  In  parallel 
with  a  voltage  monitor,  while  SMU3  acts  as  a  zero  voltage  source  In  series  with  a 
current  monitor.  Physical  connections  between  the  diode  under  test  and  the 
HP4145B  take  place  by  means  of  an  HP16058A  test  fixture.  The  fixture  and  Its 
connections  for  the  diode  case  are  depicted  In  appendix  A.  The  typical  range  of 
measurements  taken  are  shown  In  figure  3. 

Once  obtained,  measured  data  Is  submited  by  means  of  an  HPIB  (IEEE-488) 
Interface  to  a  computer  program  on  the  Hewlett  Packard  981bs  microcomputer.  The 
program  fits  a  curve  of  the  form  of  equation  2  to  the  data  yielding  values  for 
Theta  and  l,,.  The  curve  fitting  method  used  Is  described  In  appendix  B.  A 
comparison  plot  between  typical  measured  data  and  Its  resulting  fitted  curve  Is 
shown  In  figure  4. 

Bulk  Resistance 

Bulk  resistance  (RB)  Is  obtained  from  the  horizontal  difference  between  the 
curves  of  figure  4  for  high  values  of  junction  current.  Hie  migration  of 
measured  data  to  the  right  of  the  fitted  diode  equation  is  due  to  the  forward 
voltage  drop  across  RB  and  can  be  computed  from  figure  4  as, 


The  data  used  to  fit  the  curve  does  not  Include  the  last  few  points  of  measure 
data.  rhls  is  done  to  separate  the  effects  of  Rb  from  the  proper  determination 
of  parameters  Theta  and  l,,  as  well  as  legitimizing  the  calculation  of  RB  by  means 


Leakage  Resistance 


Leakage  resistance  (RS)  is  obtained  from  the  inverse  slope  of  a  curve  fit 
(app  B)  of  the  diode's  reverse  leakage  current  versus  Its  reverse  bias  voltage. 
The  following  is  the  assumed  form  to  which  the  data  is  fitted: 


The  circuit  used  to  obtain  the  reverse  bias  data  is  demonstrated  Ln  figure  5. 
Voltage  V  is  measured  across  SMU1  while  current  I  is  measured  by  SMlB.  Except 
for  the  direction  of  the  current  generated  by  SMU1,  figure  5  Is  Identical  to 
figure  2.  The  H.P.  9816S  regulates  the  measurement  taking  procedure  by  applying 
typical  currents  on  the  order  of  a  nanoarapere  and  voltages  as  high  as  40  volts. 

Diffusion  Capacitance  Constant 

Constant  KD  is  obtained  from  storage  time  measurements  collected  using  the 
circuit  of  figure  6.  A  periodic  positive  and  negative  pulse  is  applied  to  the 
diode's  junction.  Forward  and  reverse  values  of  junction  current  (Ip  and  IR, 
respectively)  are  recorded  along  with  the  time  (tg)  needed  for  the  junction  to 
discharge  and  begin  to  recover  to  its  steady  state  reverse  saturation  current 
(Ig).  The  current  seen  through  the  scope  monitored  100-ohm  resistor  is  demon¬ 
strated  ln  figure  7.  A  value  for  KD  is  obtained  using  the  following  equation 
(refs  1  and  2): 


KD  -  Theta  (- 


ln[l  +  (-r^)] 
R 


Transition  Capacitance  Equation  Constanta 


Transition  capacitance  equation  constants  Co,  PHI,  and  N  are  obtained  from  a 
curve  fit  of  transition  capacitance  (CJ)  measured  as  a  function  of  reverse  junc¬ 
tion  voltage.  The  H.P.  9816S  driven  circuit  used  to  obtain  the  desired  data  Is 
shown  in  figure  8.  The  bias  voltage  depicted  is  supplied  by  the  4145B  as 
directed  by  the  9816S.  Voltage  and  capacitance  values  as  well  as  measurement 
taking  Instructions  are  transfered  to  and  from  the  981f>S  microcomputer  by  means 
of  the  HP- I B( IEEE-488)  bus.  The  curve  fitting  procedure  discussed  in  appendix  B 
assumes  the  form  of  the  transition  capacitance  to  be 


<PHI-Vcd) 


where  V^p  represents  the  diode's  junction  voltage.  A  comparison  plot  between 
typical  values  of  measured  capacitance  and  a  curve  fit  of  the  measured  data  is 
shown  in  figure  9. 
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Diode  Node!  for  the  1N457 


The  procedures  described  in  the  previous  sections  were  implemented  in  the 
development  of  the  following  model  and  circuit  simulation: 

MODEL  DESCRIPTION 
MODEL  1N457( 1-3) 

ELEMENTS 

CD, l-2-Xl(2.3E-ll/ABS(.8-VCD)**. 874+1. 96E-5*(JIH-1.58E-10)  ) 

RB, 2- 3-3. 45 
RS, 1-2-7. 7E10 

JD, 1-2-DIODE  Q( 1 .58E-10 ,22.72) 

CIRCUIT  DESCRIPTION 

ELEMENTS 

R1  ,2-3=100 

R2  ,5-3=47 

El , l-2-TABLEl 

E2, 1-5=3 

D l, 3-6-MODEL  1N457 

RV  ,6-1  =  100 

OUTPUTS 

El  ,VRV,IRV 

CDDl ,VCDD1 

FUNCTIONS 

TABLE1 

0,4,lE-6,4,lE-6,-7,5E-6,-7 
RUN  CONTROLS 

INTEGRATION  ROUTINEpIMPLICIT 

RUN  INITIAL  CONDITIONS 

MAX  INTEGRATION  PASSES  =  1E30 

MINIMUM  STEP  SIZE  =  IE-30 

COMPUTER  TIME  LIMIT  =  1 

STOP  TIME  =  5E-6 

END 

*E0R 

8 

DBEND 

The  resulting  output  appears  In  figure  10.  The  results  closely  match  the 
measured  data  depicted  by  figure  7.  The  only  discrepency  is  the  time  needed  for 
the  diode  current  to  go  to  zero  after  its  junction  becomes  reverse  biased.  The 
difference  occurs  because  the  model  does  not  include  capacitances  which  exist 
throughout  other  parts  of  the  circuit. 


TRANSISTOR  MODEL 

The  SCEPTRE  transistor  model  (NPN)  to  be  used  is  based  on  the  Ebers-Moll 
model  (ref  1)  and  appears  in  figure  11.  The  model  parameters  are  defined  as: 


(PHI0  -  VCE 


Emitter  transition 
capacitance 


r»-  +  ^"etaE  VJE  +  IKs> 

)  E 


Emitter  diffusion 
capacitance 


- - - IT  +  Theta.  TS(JC  +  ICc) 

<PHIc  "  vcc>Nc 


Collector  transition  Collector  diffusion 
capacitance  capacitance 

CoE  ■  Constant  of  the  emitter  transition  capacitance  equation 

Cog  *  Constant  of  the  collector  transition  capacitance  equation 

PHIg  =  Emitter  base  junction  contact  potential 

PHIg  -  Collector  base  junction  contact  potential 

JE  -  Current  generator  representing  the  emitter  base  junction 

JE  *  IEg[EXP(ThetsE  VCE)  -  1] 

JC  *  Current  generator  representing  the  collector  base  junction 
JC  -  ICs[EXP(Thetac  Vcc)  -  1] 

Ng  =  Emitter  junction  grading  constant 
Nq  *  Collector  junction  grading  constant 

JN  =  Currect  generator  dependent  on  the  emitter  base  junction  current 

JN  =  01..  JE 
N 

JI  »  Current  generator  dependent  on  the  collector  base  junction  current 
JI  *  a j  JC 

IEg  *  Emitter-base  saturation  current  measured  in  the  active  region 

ICg  *  Collector-base  saturation  current  measured  in  the  active  region 

Thetag  ■  Constant  of  the  emitter  base  junction  equation 

Thetag  -  Constant  of  the  collector  base  junction  equation 

TE  *  Time  constant  of  the  emitter  diffusion  capacitance  equation 


*  t  «  ,  *  ,  Xk  If,  i  ^ **.  I'm  *  •  •"«  *' .  V/.  *  .  •  m*  *  ’  •  *•***“%  *  »  *  k  *«  *.  -  . 


6.28  Fn 

TS  =  Time  constant  of  the  collector  diffusion  capacitance  equation 


6.28  Fj 

Fn  =  The  average  ft  normal 
Fj  =  The  average  ft  inverse 

TCI  =  Charge  control  parameter  related  to  storage  time 


a  =  Forward  current  gain  =  — 

N  lE 

n  Is  entered  as  a  function  of  JE  in  the  SCEPTRE  model. 

N 

ft  =  Inverse  current  gain  =  — 

1  Lc 

n  is  entered  as  a  function  of  JC  in  the  SCEPTRE  model. 

RB  =  Base  bulk  resistance 
RC  =  Collector  bulk  resistance 
R1  =  Emitter-base  junction  leakage  resistance 
R2  =  Collector-base  junction  leakage  resistance 
Junction  Constants  and  Saturation  Currents 

Constants  Theta-  and  Theta^  represent  values  for  the  base-to-emitter  and 
base-to-collector  junction  constants,  respectively.  Similarly,  parameters  IE^ 
and  [Cg  represent  base-to-emltter  and  base-to-collector  saturation  currents. 
These  values  are  determined  In  the  same  way  as  in  the  diode  case.  The  test  setup 
used  to  obtain  the  forward  chacterlstic  of  the  base-to-emltter  junction  Is  shown 
In  figure  12,  while  the  circuit  used  to  obtain  the  characteristic  of  the  base-to- 
collector  junction  is  depicted  In  figure  13.  An  active  region  col lector-to- 
emitter  voltage  on  the  order  of  that  seen  by  the  transistor  during  application 
should  be  applied  by  SMU2  during  this  test.  In  this  test,  SMU2  acts  as  a  voltage 
source  in  series  with  a  current  monitor.  Physical  connections  between  the 
transistor  under  test  and  the  HP4145B  take  place  by  means  of  an  HP16058A  test 
fixture.  The  fixture  and  its  connections  for  the  transistor  case  are  depicted  In 
appendix  A. 


Forward  and  Inverse  Current  Gains 


Forward  current  gain  (a^)  and 
reference  1  as  follows:  ' 


inverse  current 


gain  (a  )  are 


defined  hy 


a 


N 


and 


(7  and  8) 


o  is  obtained  from  the  data  measured  by  the  circuit  of  figure  12  while  determin¬ 
ation  of  dj  is  based  on  the  results  of  the  test  setup  of  figure  13. 
o  and  are  determined  for  different  values  of  Ip  and  I,,,  respectively, 
allowing  them  to  be  described  within  the  model  in  a  tabular  way  as  functions  of 

current. 


Emitter  and  Collector  Transition  Capacitances 

As  with  the  junction  constants  and  saturation  currents  previously  discussed, 
the  emitter  and  collector  transition  capacitance  parameters  (Np,  Cop,  PHIp ,  N'^, 
Co^,  PHIq)  are  also  determined  as  in  the  diode  case.  The  test  circuits  for  each 
junction  appear  in  figures  14  and  15. 

Base  and  Collector  Bulk  Resistance 

Bulk  resistance  (RB)  is  obtained  from  the  base-to-emitter  forward  character¬ 
istic.  Its  determination  is  the  same  as  in  the  diode  case.  Collector  bulk 
resistance  will  be  assumed  to  be  typically  on  the  order  of  1  ohm.  A  more  precise 
value  for  RC  can  be  obtained  by  methods  described  in  reference  l. 

Junction  Leakage  Resistances 

The  base-to-emitter  leakage  resistance  (Rl)  and  col lector-to-base  leakage 
resistance  (R2)  are  obtained  from  the  junction  reverse  bias  characteristics.  The 
procedure  is  the  same  as  in  the  diode  case.  The  test  setups  used  are  described 
in  figures  16  and  17.  In  these  cases,  SMUl  acts  as  a  current  source  and  voltage 
meter  while  SMU3  supplies  a  voltage  of  zero  and  measures  the  junction's  current. 

Tine  Constant  of  the  Emitter  Diffusion  Capacitance  Equation 

Parameter  TE  is  the  emitter  diffusion  time  constant  and  is  introduced  in  an 
attempt  to  describe,  in  terms  of  a  time  constant,  the  capactive  effect  which  is 
observed  in  the  transistor's  emitter  junction  as  a  function  of  forward  base-to- 
emitter  voltage.  The  development  of  such  a  relationship  begins  by  expressing 
junction  voltage  in  terras  of  parameter  TE.  This  is  done  by  equating  the 
exponetlal  product  of  Theta^  and  VCp  of  the  equation  for  current  generator  JE  to 
the  familiar  exponential  terra  of  time  divided  by  some  time  constant.  In  this 
case,  the  time  constant  will  be  called  TE  and  will  produce  the  following  expres¬ 
sion  for  junction  voltage: 


1 


An  expression  for  the  emitter  diffusion  capacitance  (CEp)  as  a  function  of  TE  can 
now  be  developed.  Since  capacitance  is  defined  as 


the  derivative  of  equation  9  with  respect  to  time  can  be  subsituted  into  equation 
10  yielding  the  following  expression  for  capacitance  CEp: 


CEn  " 


=  Theta  TE  (JE  +  IE  ) 
E  S 


Theta  TE 
h 

The  above  expression  describes  the  relationship  between  capacitance  CEp  and 
parameter  TE.  Values  for  CEp  can  now  be  determined  once  parameter  TE  is 

defined.  With  this  in  mind,  an  expression  for  TE  is  next  developed  from  the 

charge  voltage  relationship  across  capacitor  CEp.  Starting  with  the  charge 
voltage  relationship  across  a  capacitor 

qE  =  CEd  Vcf^  (12) 

where  Qp  is  the  charge  on  the  junction  diffusion  capacitance  CEp  with  voltage 

Vr,  across  it.  Since  a  determination  of  TE  would  be  most  desirable  in  terms  of 

u-.p 

emitter  current,  it  is  necessary  to  rewrite  equation  12  as  follows: 

Qe  =  CEp  r.E  RE  (13) 

where  IF  is  the  emitter  current  of  the  transistor  and  RE  is  the  equivalent 
emitter  junction  resistance  seen  by  current  generator  JE  and  dimensionally  equal 
to  l/ThetaE(JE  +  rEg).  differentiating  both  sides  of  equation  13  with  respect  to 
IK  yields 


CE|>  RE  "  fhet-a”(JF-rfO  "  TB 

h  S 


The  term  —  has  the  units  of  a  time  constant  and  can  be  seen  from  equation  11  to 
dIK 

be  equal,  to  the  emitter  diffusion  time  constant  TE.  Through  the  use  of  this 
definition  and  the  equation  of  charge  continuity  (eq  13),  an  expression  for  TE 
can  he  developed  in  terms  of  emitter  current  and  time 

dQP  Q8 
3 


•-  _  A  .V  '*  % 


The  above  equation  states  that  the  base  current  of  the  transistor  is  made  up  of 
the  recombination  of  excess  carriers  every  X  seconds  and  the  time  varying  change 
in  excess  carrier  distribution.  By  separating  variables  and  then  inverting  each 
side  of  the  equation, 


lB  "  X 


u 

By  multiplying  through  by  dIR,  substituting  for  dIR  on  the  right  side  of  the 

equation,  and  then  integrating,  a  solution  for  TE  in  terms  of  iF  can  he  obtained 
as  follows: 


rBl  -5s. 
AB 


tr / 

!  djt 

I  dQ~c 


Limit  t  represents  the  time  it  takes  for  the  emitter  current  to  rise  form  27%  to 
73%  of  its  final  value.  Ifl  is  assumed  to  go  from  zero  to  Igl  in  a  stepwise 

Qb 

manner.  Since  - -  represents  the  d.c.  component  of  the  base  current,  it  can  be 

XB 

subsituted  by  the  d.c.  emitter  current  divided  by  $  +  1  as  follows: 


(B  +  1) 


(8  *  1)  I„  -  I„ 


27%I_y  6 

E 

B  is  the  average  B  for  the  27%  to  73%  collector  current  transition  experienced. 

Values  of  B  for  different  values  of  collector  current  can  be  calculated  from 

the  a  results  obtained  from  equation  7  and  the  following  expression  for  g  . 

N 


(i  -«„> 


Integrating,  the  following  is  obtained: 


Tjjf  =  -(P  +  1)  ln[(B  +  1)  IB1  -  IEJ 


;  27%It 


By  evaluating  equation  20  over  the  limits  indicated  and  then  substituting  IF 
for  (B  +  D  Igl*  t^le  f°How^ng  expression  for  TE  can  be  developed: 


TE 


(21) 
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(3  +  0  In 


0.73 

0.27 


t 
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3  +  1 


Therefore,  the  value  of  TE  for  different  saturation  currents  can  be  obtained  by 
measuring  the  time  necessary  for  the  collector  current  or  collector-to-emitter 
voltage  of  the  transistor  to  go  from  27%  to  73%  of  its  final  value.  Therefore, 
TE  is  expressed  within  the  model  in  a  tabular  manner  as  a  function  of  current 
generator  JE  and  is  obtained  from  the  test  setup  of  figure  18.  An  example  of  the 
oscilloscope  trace  generated  by  the  test  circuit  of  figure  18  is  shown  in  figure 
19. 


Tine  Constant  of  the  Collector  Diffusion  Capacitance  Equation 

The  time  constant  (TS)  of  the  collector  diffusion  capacitance  equation  is 
determined  by  observng  the  depletion  of  excess  charge  in  the  transistor's  base 
region  as  it  is  brought  out  of  saturation.  By  using  the  equation  of  charge  con¬ 
tinuity,  reference  3  develops  an  expression  for  the  transistor's  storage  time 
constant  as  a  function  of  base  current  and  collector  saturation  current.  The 
expression  is  as  follows: 


dQ 


BS 


dt. 


,(IBl  IB2). 


(22) 


In  the  above  expression,  IB)  represents  the  base  current  while  the  transistor  is 
in  saturation;  lR.,,  the  base  current  used  to  bring  the  transistor  out  of  satura¬ 
tion;  and  Tr,  the  collector  saturation  current.  Parameter  t  is  the  storage  time 
needed  to  deplete  the  base  of  excess  charge  carriers  and  begin  to  bring  the 
transistor  out  of  saturation.  Values  for  I^j  ,  ,  and  ts  can  be  obtained  by 

using  t he  circuit  of  figure  18.  The  method  used  to  obtain  a  measurement  for  t 
is  shown  in  figure  19.  Parameter  t,  is  measured  from  the  time  when  the  transis¬ 
tor's  base  is  at  its  minimum  to  the  time  its  collector  voltage  reaches  10%  of  its 
final  value.  The  10%  point  of  reference  is  used  to  make  it  easier  to  obtain  a 
reading  while  introducing  only  a  negligible  amount  of  error  into  the 
measurement.  The  result  of  equation  22  and  the  inverse  mode  model  components 
IC,  CC,  and  .11  of  figure  10  are  used  to  explain  saturation  region  behaviour  and 
develop  time  constant  TS  in  terms  of  current  generator  JC.  TS  is  defined  as: 


TS  = 


dO 
d  f 


BS 


/  %  _1«2_\ 
VB  -  h- 


(23) 


Parameters  3^,  8  ^  ,  T^,,  and  o  are  all  functions  of  the  operating  point  and  are 

obtained  by  a  SCEPTRE  simulation  (ref  1).  The  simulation  consists  primarily  of  a 
constant  capacitance  steady  state  model  of  the  transistor  which  uses  all  of  the 
parameters  obtained  up  to  this  point.  The  model  Is  excited  in  a  common  emitter 
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circuit  at  different  values  of  collector  current.  The  results  produce  values 
from  which  parameter  TS  can  he  expressed  In  tabular  form  as  a  fun  t  ion  it  l'i  • 
model  current  generator  JC. 


k  Transistor  Model  for  the  2N2222 

The  procedures  described  In  the  previous  sections  were  implemented  i  1  t  ue 
development  of  a  model  for  the  2N2222  transistor.  Tabulated  test  results  used  n; 
the  determinat Ion  of  parameters  TE  and  TS  appear  li  tables  \  and 
respectively.  A  listing  of  the  model  and  the  SCEPTRE  input  used  to  simulate  toe 
circuit  of  figure  IB  follow: 

MODEL  DESCRIPTION 
MODEL  2N2222( B-E-C) 

ELEMENTS 

CE , 1 -E-Q 1 ( 4 . 5E- l I ,.75,VCK,.78,28,T3( JE) ,JK,2.2K- I  I) 

CC  ,  1-2-Q l ( 2 .9E-1  1  ,.75  ,  VCC  ,1.03,28.4,T4(  IC)  ,  IC.2.8E-1  1 

RB. B-1-47 

RC. C-2-.5 

R1 .1-E-3.9E10 
R2.1-2-1.1E10 

JE, 1-E-DIODE  Q<  2 . 2E- 1 1 ,28) 

JC  ,  1-2-DIODE  Q(2.8E-1 1 ,28.4) 

JN, 2-1 -TABLE l ( JE)*JE 
J l , E- 1 -TABLE2( JC) * JC 
FUNCTIONS 

Q1(A,B,C,D,E,F,G,H)-(A/ABS( B-C)**D+E*F*(G+H) 

TABLEl-0,.9925,4. IE-4,.9925,1. IE-3,.9928 .2.35K-3, .9928  ,  l  ,:i4K-  ’ , 

.9926 , l . I 5E-2 , . 9913, 1. 37 E- 2, .992, 2 E- 2, .992 

TABLE2-0 , 5 .036E-2  ,2 .0 1 9E-6 ,5 .<  36E-2  ,2 .6  7  E-  S  , .  102  8  ,  S  .  HE-  5  ,  .  H8n  , 

8.  1  36E-  5  , .  I  66 , 1  .  2 44E- 4 , .  I  9  7 , 2 . 6  E- 4  , .  2  6  I  5 , 3 E-  4 , .  2  8  l  5 

TABLE3-0  , .  5E-9 ,1  5E- 3  , .  5F.-9  ,2'>E-  3 ,  .4  2  K-9 ,30F-  J ,  .SK-9  ,4uK-  l,  .8  Ur  -‘i, 

1 , .833E-9 

TABLF.4-0 ,105.  7E-9  ,5.  34E-4 , 1  OS.  7E-9 , 1  .1  I  7E-  J.l  I  7K-R  .2 . 44K-  i,>8.ot-9 
1 .58.6E-9 

CIRCUIT  DESCRIPTION 
ELEMENTS 
Rl .2-3-1E3 
R2  ,5-4-100 
El  , t -2-TABLE  1 
E2  ,1-5-4 

Tl ,3-1-4-MODEL  2N2222 
J1  ,4-1=0 
OUTPUTS 
V 11 ,E  l 


FUNCTIONS 
TABLE l 

0,-1. 6 .1200E-9  ,-1.6 , 1  280K-9,  I  .6  .2480E-9 ,1  . 6 , 2  8  JOE- 9  ,  -  1 .  h  ,2  8;  iOF-  >» ,  -  l  .  h 
RUN  CONTROLS 

I  INTEGRATION  ROUT  I NE-  l  MP  L I C  IT 

RUN  INITIAL  CONDITIONS 

MAX  INTEGRATION  PASSES  -  IE10 


1  I 


MINIHUM  STEP  SIZE  -  IE- 30 

CO*>UTER  TIME  LIMIT  -  1 

STOP  TIME  -4000E-9 

END 

*EOR 

8 

DBEND 

Its  resulting  output  appears  in  figure  20.  The  results  closely  match  the 
measured  data  depicted  by  figure  19. 


CONCLUSIONS 

The  results  obtained  in  this  report  confirm  the  modeling  procedures  used  by 
Howers  and  Sedore  and  demonstrate  a  method  for  successfully  automating  those 
procedures  using  the  Hewlett  Packard  98lhS  microcomputer  and  4 1 4  SB  semiconductor 
parameter  »na  l  yr.er. 


Table  1.  Transistor  current  rise  time  parameters  for  the  2N2222 
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igure  4.  f'urve  fir  of  the  forward  characteristic  of  a  1N457 
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Reverse  diode  characteristic  test  circuit 
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Figure  /.  Diode  (1N457)  reverse  recovery  waveform  from  the 

test  setup  of  the  diffusion  capacitance  test  circuit 


Figure  8.  Transition  capacitance  test  circuit 
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Simulated  diode  current  of  the  diffusion 
capacitance  test  circuit 
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Physical  connections  between  the  devices  under  test  and  the  HPAlAbB  semicon¬ 
ductor  parameter  analyzer  are  discussed  In  this  appendix.  Connections  are  made 
through  the  use  of  an  HP16058A  test  fixture.  The  wiring  layouts  used  in  the 
diode  and  transistor  tests  are  demonstrated  In  figures  A-l  and  A-2.  A  general 
list  of  the  equipment  used  In  this  report  Is  shown  In  table  A-l. 

While  acting  as  an  Interface  between  the  HP414SB  and  the  device  under  test 
(DUT),  the  test  fixture  also  provides  a  means  of  easily  modifying  circuit  connec¬ 
tions  with  the  use  of  a  connection  switch.  The  switch  in  figure's  A-l  and  A-2 
connects  the  center  column  nodes  (COMM)  to  the  nodes  of  either  columns  1  or  2. 

When  performing  diode  related  tests,  the  layout  of  figure  A-l  is  used.  The 
connection  switch  is  moved  to  position  1  when  the  diode’s  forward  characteristic 
and  leakage  resistance  are  being  determined.  The  switch  Is  moved  into  position  > 
when  transition  capacitance  tests  are  being  performed,  hut  the  diode  must  he 
removed  from  the  te«t  fixture  in  this  case  and  connected  directly  to  the  I. OR 
meter. 

Transistor  tests  are  made  using  the  layout  of  figure  A-2.  Once  again,  the 
connection  switch  is  used  to  modify  test  procedures  without  disturbing  the  test 
fixture.  The  switch  is  placed  In  position  1  when  performing  the  forward  NORMAL 
mode  characterist  ic  test  and  in  position  2  when  determining  the  INVKKSK  mode 
character  1st Ic.  Transition  capacitance  tests  are  made  using  the  layout  of  figure 
A-l  while  remembering  that  the  diode  depleted  represents  the  translst>r  Junction 
under  test.  Junction  leakage  resistance  tests  use  the  layout  of  figure  A-2  with 
modifications  stipulated  by  base-to-emi t te r  and  base-to-col lector  reverse 
character  1st Ic  test  setups  described  earlier  in  this  report. 


\  < 


td 


Table  A-l.  Equipment  list 

1.  Hewlett  Packard  9816S  microcomputer. 

2.  Hewlett  Packard  4145B  sealcoductor  parameter  analyzer. 

3.  Hewlett  Packard  4262A  LCR  meter. 

4.  Lambda  LE103  PM  power  supply. 

5.  Tektronix  5441  storage  oscilloscope. 

6.  Tektronix  FG504  function  generator. 
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Figure  A- 1  .  Test  fixture  (H.P.  lhOSHA)  di 


APPENDIX  B 

LEAST-SQUARES  CURVE  FITTING 
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The  curve  fitting  method  used  in  this  report  is  known  as  the  least-squares 
method.*  It  is  based  upon  the  concept  that  the  best  line  fitted  to  a  set  of  data 
points  is  the  one  for  which  the  sum  of  the  squares  of  the  deviations  between  the 
points  and  the  line  is  a  minimum. 

The  deviations  are  squared  to  eliminate  the  possible  cancelling  of  positive 
and  negative  deviations  at  different  points  which  could  give  a  false  impression 
of  a  good  fit. 

This  report  maps  all  of  the  expressions  to  be  used  for  fitting  purposes  to 
the  form  of  the  straight  line,  Y  =  A  +  BX.  The  sum  of  the  squares  of  the  devia¬ 
tions  may  be  expressed  as 


Z(Vyo: 


(A+Bx  -  y  )“ 


(B-L) 


where  yc  *  A  +  BxQ  is  the  calculated  ordinate  for  the  value  of  x  =  xQ. 

The  values  of  the  constants,  A  and  B,  for  which  this  expression  is  a  minimum 
may  be  determined  by  finding  the  partial  differentiation  of  equation  B-l  with 
respect  to  each  constant  and  then  setting  it  equal  to  zero. 


For  A: 


For  B: 


(A  +  Bx  -  y  )  *  0 

o  o 


2  (A+Bx  -  y  )  =  0 
/  .  o  J  o 


(A+Bx  -  y  )  =  0 
o  J  o 


T~ '  ( A  +  Bx  -  y)2  =  0 
3  B  O  O 


2  (Ax  +  Bx  -  x  y  )  =  0 
L  ,  o  o  o  o 


(Ax  +Bx  -  x  y  )  =  0 
o  o  oJ  o 


These  two  equations  can  be  solved  simultaneously  for  unique  values  of  the 
quantities  A  and  B  since  values  of  xQ  and  yQ  are  available  from  the  original 
data.  The  two  equations  may  be  written  in  the  form: 


*  Barnas,  Robert  E. ,  Cuvit  II  -  A  Curve  Fitting  Program  with  Plotting  Options, 
Management  Information  Systems  Directorate,  Picatinny  Arsenal,  Dover,  NJ. 
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x0y0  ■  Axo +  L*\ 


Factoring  A  and  B  from  the  terms  gives 


*  N  A  +  B  '<•  x. 

J  a  r 


y.x.  =  A  x  +  B  >  x 
/  o  o  u 0  /.  .  o 


where  N  =  number  of  data  points  and  finally 
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The  manner  in  which  assumed  data  forms  were  mapped  into  the  straight  line  form  is 
described  in  table  B-l. 
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Table  B-l.  Straight  line  equation  mappings  used  in  the 
least-squares  fit  of  semiconductor  data 


Equation  no. 

Original  form 

Desired  form 

Mapping 

2 

I=ls  Exp(Theta  V^d) 

V 

Y=A+BX 

Y=ln( I) 

A=ln(Is) 

B=Theta 

X=V 

4 

I=  RS 

Y=A+BX 

Y=I 

A=0 

B=  RS 

X=V 

6 

CJ=  Co 

N 

((PHI-VJ)  ) 

Y=A+BX 

Y=ln(CJ) 

A=ln(Co) 

B=-N 

X=ln( PHI-VJ) 
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